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The detection and characterization of paramagnetic species by electron-spin 
resonance (ESR) spectroscopy is widely used throughout chemistry, biology, and 
materials science [1], from in-vivo imaging [2] to distance measurements in spin- 
labeled proteins [3]. ESR typically relies on the inductive detection of microwave 
signals emitted by the spins into a coupled microwave resonator during their 
Larmor precession — however, such signals can be very small, prohibiting the 
application of ESR at the nanoscale, for example, at the single-cell level or on 
individual nanoparticles. In this work, using a Josephson parametric microwave 
amplifier combined with high-quality factor superconducting micro-resonators 
cooled at millikelvin temperatures, we improve the state-of-the-art sensitivity of 
inductive ESR detection by nearly 4 orders of magnitude. We demonstrate the 
detection of 1700 bismuth donor spins in silicon within a single Hahn [4] echo 
with unit signal-to-noise (SNR) ratio, reduced to just 150 spins by averaging a 
single Carr-Purcell-Meiboom-Gill sequence [5]. This unprecedented sensitivity 
reaches the limit set by quantum fiuctuations of the electromagnetic field instead 
of thermal or technical noise, which constitutes a novel regime for magnetic 
resonance. The detection volume of our resonator is ~0.02 nl, and our approach 
can be readily scaled down further to improve sensitivity, providing a new and 
versatile toolbox for ESR at the nanoscale. 

A wide variety of techniques are being actively explore to push the limits of sensitivity of 
ESR to the nanoscale, including approaches based on optical [6, 7] or electrical [8, 9] detec¬ 
tion, as well as scanning probe methods [10, 11]. Our focus in this work is to maximise the 
sensitivity of inductively detected pulsed ESR, in order to maintain the broad applicability 
to different spin species as well as fast high-bandwidth detection. Pulsed ESR spectroscopy 
proceeds by probing a sample coupled to a microwave resonator of frequency uq and quality 
factor Q with sequences of microwave pulses that perform successive spin rotations, trigger¬ 
ing the emission of a microwave signal called a spin-echo whose amplitude and shape contain 
the desired information about the number and properties of paramagnetic species. The spec¬ 
trometer sensitivity is conveniently quantified by the minimal number of spins Nmin that 
can be detected within a single echo [4]. Conventional ESR spectrometers use 3D resonators 
with moderate quality factors in which the spins are only weakly coupled to the microwave 
photons and thus obtain a sensitivity of N^in 10^^ spins at T = 300 K and X-band fre- 
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quencies {u:q/2t[ ~ 9 — 10 GHz). To increase the sensitivity, micro-fabricated metallic planar 
resonators with smaller mode volumes have been used, resulting in larger spin-microwave 
coupling [12, 13]. Combined with operation at T = 4K and the use of low-noise cryogenic 
amplifiers and superconducting high-Q thin-hlm resonators, sensitivities up to Nmin ~ 10^ 
spins have been reported, which represents the current state-of-the-art [14-16]. 

Further improvements in the sensitivity of ESR spectroscopy can be obtained by cooling 
the sample and resonator down to mK temperatures that satisfy T Uuo/kB at X-band 
frequencies. As a result, both the spins and the microwave held reach their quantum ground 
state, which is the optimal situation for magnetic resonance since the spins are then fully 
polarized and thermal noise suppressed. The noise in the emitted echo signal is essentially 
due to vacuum quantum huctuations of the microwave held, with a dimensionless spectral 
power density Ueq = S{u)/{hu) = 1/2, possibly supplemented by extra noise n* due to 
spontaneous emission of the spins (see Supplementary Information). The total noise spec¬ 
tral density in the detected signal n = Ueq + -f Uamp however also includes the added noise 
namp of the hrst ampliher of the detection chain; benehting from the low noise ahorded by 
low temperature operation thus requires nearly noiseless amplihers at microwave frequen¬ 
cies, as were recently developed in the context of superconducting quantum circuits. These 
Josephson Parametric Amplihers (JPAs) are operated at mK temperatures, have a band¬ 
width of up to ~ 100 MHz, and a low saturation input power (typically 1 — 10 fW) [17, 18]. 
They have been shown to add the minimum amount of noise permitted by quantum me¬ 
chanics [19]: Uamp = 0.5 when both held quadratures are equally amplihed (non-degenerate 
mode) [17], and Uamp = 0 when only one quadrature is amplihed (degenerate mode) [20]. 
JPAs have been used so far for reading-out the state of superconducting qubits [21], the 
motion of nanomechanical oscillators [22] and the charge state of a quantum dot [23], as 
well as for high-sensitivity magnetometry [24]. Here we show that they are also well suited 
to amplify the weak and narrow-band signals emitted by small numbers of spins, with the 
ultimate sensitivity allowed by quantum mechanics, enabling us to demonstrate a 4 orders 
of magnitude improvement in sensitivity over the state-of-the-art. 

We use an ensemble of bismuth donors implanted over a 150 nm depth into an isotopi- 
cally enriched silicon-28 crystal, on top of which we pattern a superconducting aluminium 
thin-hlm micro-resonator consisting of an interdigitated capacitor in parallel with a wire 
inductance (see Fig. 1 for a sketch of the setup). Due to this geometry, the microwave held 
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FIG. 1. Experimental setup and spin system, (a) The aluminium microwave resonator with 
frequency wq consists of an interdigitated capacitor in parallel with a 5/rm-wide wire inductor, 
fabricated on a Bi-doped ^®Si epi-layer. (b) The sample is mounted in a copper box, thermally 
anchored at 12 mK, and probed by microwave pulses via asymmetric antennas coupled to the 
resonator with rate ki = 1.2 • 10^s“^ and K 2 = 5.6 • 10^s“^. A magnetic field Bq is applied parallel 
to the resonator inductance. Microwave pulses at ujq are sent by antenna 1, and the microwave 
signal leaving via antenna 2 is directed to the input of a Josephson Parametric Amplifier (JPA). 
The JPA is powered by a pump signal at iOp fs 2a;o, and its output is further amplified at 4K 
by a High Electron-Mobility Transistor amplifier, followed by amplification and demodulation at 
room-temperature, yielding the two field quadratures (c) Energy levels of Bi donors 

in Si, expressed in frequency units (see spin Hamiltonian in the Supplementary Material), (d) 
ESR-allowed transitions in the low-field limit. Eor Bq < 8mT, the \F,mF) = |4, —4) —>■ |5,—5) 
and |4, —3) —)• |5, —4) transitions cross the resonator frequency at respectively Bq = 5 and 7mT. 
(e) Measured resonator transmission coefficient |S' 2 i| (red circles), yielding t(Jo/2vr = 7.24GHz and 
a total quality factor Q = 3 - 10^ (red curve is a fit), (f) The JPA can be characterized via a direct 
line bypassing the resonator, yielding a gain, in non-degenerate mode, of G > 20 dB above a 3 MHz 
bandwidth. Circles are experimental data, curve is a Lorentzian fit. 
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Bi cos oj^t couples only to the Nsi — 4 ■ 10^ implanted Bi atoms located in the area below 
the wire. The sample is inserted inside a copper box to suppress the resonator radiative 
losses while enabling to probe its transmission by capacitive coupling to input and output 
antennas. In this well-controlled environment the resonator reaches a loaded quality factor 
Q = 3 X 10^ for a frequency ujq/2ti = 7.24GHz (see Fig. le). Microwave pulses at uq are 
applied to the cavity input; the output signal (including the echoes emitted by the spins) 
is directed towards the input of a JPA [18] with power gain G up to ~ 23 dB at Uq when 
powered by a pump microwave signal at a frequency Up ~ 2uo, ojp = 2ujq corresponding to 
the degenerate mode of operation [18]. The JPA output is then further amplihed by a semi¬ 
conducting HEMT ampliher at 4K, and hnally demodulated at frequency ujq yielding time 
traces for the two quadratures I(t),Q(t) (see Supplementary Information for more details). 
The sample and JPA are cooled at 12 mK in a dilution refrigerator. 

An in-plane magnetic held Bq is applied parallel to the sample surface along the resonator 
inductance. Neutral bismuth donors in silicon have a S' = 1/2 electron spin coupled by a 
strong hyperhne interaction term A~^ ■ ”/^ to the / = 9/2 nuclear spin of ^°®Bi [25, 26], 
with A/h = 1.48 GHz. In the low-held regime, the 20 electro-nuclear energy states are best 
described by their total angular momentum ^ and its projection mp — they 

can be grouped in a F = 4 ground and a F = 5 excited multiplet separated by a frequency 
of 3A/h = 7.38GHz in zero held (see Fig. 1). With the chosen orientation of Bq, the Bi 
microwave held generated by the resonator is perpendicular to the spin quantization axis 
and only transitions obeying \AmF\ = 1 have a signihcant matrix element (see Fig. Id) for 
Bq < 10 mT. Their frequency in the ~ 7.3 — 7.5 GHz range makes Bi:Si an ideal system 
for coupling to superconducting aluminum resonators which can withstand only helds below 
~ 10 mT. 

In our case, the |F, 771 ^) = |4, —4) —)■ |5, —5) and |4, —3) —)■ |5, —4) transitions are 
expected to be resonant with uq at Fq = 5 and 7 mT respectively; corresponding peaks 
in the integrated spin-echo signal (of duration Tp ^ 20p.s) are indeed measured as shown 
in Fig. 2a-c. Each transition consists of two sub-peaks, with an inhomogeneous linewidth 
F/27r = 2 MHz. We attribute this sub-structure to the diherential strain [27] acting on the 
Bi atoms lying just under the wire versus those around it (see Supplementary Information). 
We will focus in the following on the |4, —4) —)■ |5, —5) transition for the spins lying under 
the wire, at Bq = 5.18 mT. Well-dehned Rabi oscillations are observed in the integrated echo 
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FIG. 2. Sample characterization, (a) Hahn-echo sequence (top), triggering the emission of 
an echo (bottom). Plotted are the demodulated quadratures I{t) (green squares) and Q{t) (red 
diamond), as well as the echo amplitude A{t) = + <3(t)^ (blue circles), from which the echo 

quadrature area Xe = (with X = I, Q) and amplitude area Ae = A{t)dt 

are extracted. The data were taken for Bq = 5.2 mT. (b) Normalised amplitude echo area as a 
function of the refocusing pulse amplitude A^^ (rescaled by the amplitude needed for a vr pulse) 
showing Rabi oscillations. Blue circles are data points, red curve is an exponentially damped sine 
fit. (c) Amplitude echo area (blue circles joined by dashed lines) as a function of magnetic held 
Bq showing two principal resonances, each split into a doublet due to the effect of strain on the 
donors below and next to the aluminium wire inductor, (d) As the total time 2r between the 
initial tt/2 pulse and the echo is increased, the recovered Q quadrature echo area decays with an 
exponential behaviour (red curve is a ht), yielding a spin coherence time T 2 = 8.9ms. (e) The 
inversion recovery sequence (see inset) is used to measure the spin relaxation time Ti = 0.37 s. Red 
curve is an exponential ht to the experimental &ata (blue circles). 













signal as a function of the refocusing pulse amplitude (see Fig. 2b), with a 100kHz Rabi 
frequency for a remarkably low input power of 3pW [15]. The decay of the integrated echo 
signal as a function of the total delay 2r between the initial 7r/2 pulse and the echo is well 
htted by an exponential decay with a time constant T 2 = 10 ms, a typical coherence time 
for Bi Si [28] (see Fig. 2d). The energy relaxation time Ti is measured by the inversion 
recovery method to be Ti = 0.3s (see Fig. 2e), allowing us to use a IHz repetition rate 
throughout this work. 

The spectrometer sensitivity is estimated by measuring the SNR of a single echo. The 
JPA is operated in the degenerate mode, with the phase of the pump signal chosen such 
that the echo signal is entirely on the amplihed quadrature. With these optimal settings, 
the amplitude SNR of the echo shown in Fig. 3a is found to be 7 ± 1, one order of magnitude 
larger than the SNR obtained in the same conditions but with the JPA pump turned off so 
that it simply reflects the echo signal. This improvement is consistent with a noise reduction 
from n ~ 50 (with JPA off) down to n ~ 0.5, thus close to the quantum limit, and with 
calibration measurements performed on the JPA itself (see Supplementary Information). 

Of all the neutral Bi donors within the resonator mode volume, only those whose fre¬ 
quency lies within the resonator linewidth k = ojq/Q and which are in the |4, —4) state 
contribute to the echo signal. A rough estimate of the number of spins is therefore obtained 
as A^Bj(K/P)/9 = 4-10^, an overestimate given that only a fraction of implanted atoms 
shows a magnetic resonance signal due to either crystal damage or to donor ionization [29]. 
For a more accurate determination, the time-dependent absorption of a microwave pulse 
at Uq recorded and htted to a simple model (see Fig. 3b and Supplementary Information) 
allows us to obtain an abolute calibration of the spin density. A whole spin-echo sequence 
is then measured and simulated (see Fig. 3c); the quantitative agreement with the observed 
echo amplitude establishes (from the simulations) that 1.2 x 10"^ spins are excited during the 
sequence. This implies a ~ 30% yield between number of implanted atoms and of neutral 
donors, compatible with previous reports [29]. 

Overall, the spectrometer can therefore detect down to Nmin = 1.2 x 10^/7 = 1.7 ■ 10^ spins 
with a signal-to-noise of unity in a single Hahn-echo, and has a corresponding sensitivity of 
1.7 ■ 10^ spins/v^Hz given the 1 Hz repetition rate. This 4 orders of magnitude improvement 
over the state-of-the-art is in qualitative agreement with the prediction of a simplihed model 
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FIG. 3. Spectrometer sensitivity, (a) Echo signal I{t) without JPA (red curve) and with JPA 
in degenerate mode and 23dB gain (blue curve), averaged 10 times. The red curve was rescaled 
by the amplifier amplitude gain for comparison with the JPA-on curve. The JPA-on (JPA off) 
echo has a signal-to-noise ratio of 22 ± 3 (2 ± 0.5), which translates into a single-echo signal-to- 
noise of 7 ± 1 (0.6 ± 0.15). (b) Time dependent amplitude of a 500/rs pulse at wq, averaged 1000 
times, showing Rabi oscillations (circles). A simulation (curve) is used to estimate the number of 
spins contributing to the absorption, (c) Measured microwave amplitude (circles) during an entire 
Hahn echo sequence, with the JPA turned off in order to avoid any saturation effect. A simulation 
(curve) uses only the number of spins extracted from (b) and shows quantitative agreement with 
the measurements. These simulations indicates that the '7r/2 pulse acts on 1.2-10^ spins. 

(see Supplementary Information) ~ ^ being the coupling constant of a single 

spin to the resonator microwave held, estimated for our geometry to be 5'/27r = 55 Hz, which 
yields = 400 spins. The sensitivity can be further improved with a Carr-Purcell- 

Meiboom-Gill pulse sequence, adding m vry pulses after the hrst echo in order to recover 
m echoes instead of a single one, yielding an increase in signal-to-noise of ~ [5]. The 

applicability of this technique depends on factors such as the spin coherence time T 2 of the 
sample and the echo duration Te — for our ^®Si;Bi sample up to 600 echoes are obtained, 
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FIG. 4. Further sensitivity improvement with the Carr-Pnrcell-Meiboom-Gill (CPMG) 
pulse sequence, (a) A spin echo generated by any pulsed ESR experiment can be refocused by 
a train of vr pulses with rotations axes oriented along the echo phase direction, and thus used 
to enhance the SNR in a single shot, (b) Time-dependence of the measured echo amplitude as 
a function of the echo number, m, over the course of a single sequence comprising 650 vr pulses 
separated by r = 200/is. (c) Three time traces (circles) of echoes number 1, 100, 600 are explicitly 
shown from (b). Solid line shows the average over all 650 echoes, (d) Decay of the echo area 
as a function of the total delay between the initial 7r/2 pulse and the echo (circles), fitted with an 
exponential decay (curve) of time constant Tcpmg = 71 ms. (e) SNR improvement as a function 
of the number m of vr pulses within the CPMG sequence (circles), showing a tenfold improvement. 
In the absence of decoherence, the SNR should follow ^/m (green curve); with decoherence (red 
curve) the SNR levels off, and eventually decays for higher m. 
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as shown in Fig. 4, with a corresponding tenfold increase of the SNR and an unprecedented 
sensitivity of 150 spins in a single shot or 150 spins/V^Hz. 

A wide range of species including molecular magnets, Gd spin-labels and high-spin defects 
in solids, can be studied by ESR at low magnetic helds using the A1 thin-hlm resonator 
demonstrated here. Operation in larger magnetic helds (~ 0.3 T) would enable the most 
general application of this method to other spin species and could be achieved by fabricating 
the micro-resonator from higher critical held superconductors such as Nb [15] or NbTiN [30]. 
Our results thus open the way to performing ESR spectroscopy on nanoscale samples such 
as single cells, small molecular ensembles, nanoparticles and nano-devices. We predict a 
further 2 orders of magnitude sensitivity enhancement is possible by reducing the resonator 
transverse dimensions down to the nanometric scale, which would then be sufficient for 
detecting individual electron spins. 
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Supplementary Material: Reaching the quantum limit of 
sensitivity in electron spin resonance 

I. EXPERIMENTAL DETAILS 

Bismuth implanted sample 

The sample consists of a natural silicon (100) substrate on which a 700nm-thick isotopi- 
cally enriched 99.95% ^®Si epitaxial layer was grown. Bismuth dopants were subsequently 
implanted into the epitaxial layer and activated by thermal annealing (see [29] for more 
details). The implantation prohle, measured via Secondary Ion Mass Spectroscopy (SIMS) 
is shown in SI of the main text. The activation step consists in an anneal to SOO®^ for 
20 min under nitrogen atmosphere. An electrical activation yield of 60% has been measured 
using a Hall effect measurement system under similar implantation conditions [29]. 
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FIG. SI. Bismuth implantation profile, measured via Secondary Ion Mass Spectroscopy (SIMS) 

A 50-nm-thick aluminum resonator was deposited on top of the sample using a standard 
lift-off process. The resonator is a lumped element constituted by a 5pm-wide wire and an 
interdigitated capacitance of 12 50/im-wide hngers spaced by 50 pm. The wire is 730 pm 
long yielding a resonance ujq/2t: = 7.26 GHz. 

Measurement setup 

The detailed microwave setup is shown in supplementary Fig. S2. The sample is enclosed 
in a box made of oxygen-free-high-conductivity-copper, whose lowest resonance mode is at 
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8.4 GHz. Its role is to suppress the resonator radiative losses, while enabling its coupling to 
the input and output antennas with rates Ki (input) and K 2 (output). The values of Ki and 
K 2 as well as the resonator internal losses kl are determined experimentally by measuring 
the complete resonator scattering matrix with the spins far from resonance, and htting this 
matrix to the known resonator input-output formulas (see [31] for instance). The results 
are shown in Table I. As can be seen in Table I, the experiment is in the so-called critical 
coupling regime where the internal losses kl are approximately equal to the external coupling 
Ki + K 2 . The asymmetry between ki and H 2 was purposedly chosen to be large (~ 5) so as 
to ensure that the majority of the photons emitted by the spins would be collected by the 

output antenna. The resonator total quality factor is Q = --= 3 x 10®. 

+ K2 + Kl 


UJ'p 

2^ 

Q 

Kl 

K2 

Kl 

7.24 GHz 

3 X 10® 

2.1 X 10®s-i 

9.2 X 10®s-i 

12 X 10®s-i 


TABLE I. 

During a spin-echo sequence, a series of microwave pulses at frequency cuq is sent on the 
input line of the cavity (ki, green wire on the hgure). The transmitted signal is then routed 
to the Josephson Parametric Ampliher (JPA) via a circulator and is further amplihed hrst 
by a low-noise HEMT ampliher at the 4K stage then at room temperature, before being 
hnally demodulated by mixing with the local oscillator at frequency ojq. 

The JPA can be studied and tuned independently from the cavity by an additional 
input line (brown) coupled via a 20dB coupler to the cavity output line. Its design and 
operation have been described in detail in [18]. It consists of a lumped element resonator 
formed by an interdigitated capacitance, a geometrical inductance and an array series of 8 
superconducting quantum interference devices (SQUIDs). The SQUID array acts as a hux- 
tunable inductor that allows the resonator frequency to be tuned over a 400 MHz frequency 
range by passing a dc current through an on-chip antenna. The ampliher is parametrically 
pumped by modulating the hux threading the SQUIDs at a frequency Up close to twice 
the resonator frequency. The JPA can be operated either in phase sensitive mode where 
Up/2 = uq or in non-degenerate mode. In the latter mode, around 7.3 GHz, a non-degenerate 
gain of 23 dB can be obtained with the appropriate pump power as shown in Fig. 1 of 
the main text. Saturation of the JPA occurs for a typical input power of —130 dBm. In 
all the data presented in this work, the echo signals emitted by the spins are below the 
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FIG. S2. Measurement setup. 

amplifier saturation threshold. The detuning between the pump and the signal is chosen to 
be ~ 500 kHz, and the demolutated signal is hltered at 100 kHz to suppress the idler. In 
phase-sensitive mode, only one signal quadrature is amplihed, and the JPA has an additional 
6dB gain. The quadrature is chosen by tuning the relative phase of the pump and signal 
sources. In addition, the JPA pump is pulsed via the microwave source internal switch (not 
shown on schematic= so as to generate gain only during the emission of an echo signal. This 
is done to reduce the effective pump power brought to the mixing chamber plate and thus 
bring the cryostat temperature from 20 mK with a continuous pump signal to 12 mK. 

A double circulator is used to prevent interferences between the cavity and the JPA. 
Another double circulator is needed at the JPA output to ensure both the routing of the 
signal and the isolation from the thermal photons travelling down the output line (red line). 
Leakage of the 14.5 GHz pump to the resonator and the spins is suppressed by a 4 — 8 GHz 
bandpass hlter inserted between the cavity and the JPA. Each input line is attenuated by 
20dB at 4K and 20dB at 12mK to thermalize the electromagnetic held, and hltered by low- 
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pass filters containing infra-red absorptive material in order to minimize losses due to out- 
of-equilibrium quasi-particles generated in the superconducting thin-film. Both the cavity 
and the JPA are magnetically shielded. The JPA is enclosed in an 3-cm-wide aluminium 
box surrounded by a 1-mm-thick cryoperm material, the whole being placed inside a 20- 
cm-long yU-metal cylindrical shield. The magnetic held Bq is applied parallel to the sample 
surface with an arbitrary angle 9 with respect to the resonator axis, by using two orthogonal 
Helmholtz coil pairs that can provide up to 10 mT and have been calibrated in a previous 
experiment. 

Pulses are shaped by a microwave switch in series with the microwave source internal 
gate. The relative phases of the pulses are controlled by analog phase modulation. Every 
control signal is generated by an arbitrary wavefrom generator (AWG). In order to suppress 
any offset in the detection chain without time-consuming calibration, every pulse sequence 
is repeated twice with opposite phases on the 7r/2 pulses. This phase cycling protocol yields 
two echo signals with opposite phases taken in the same conditions : the offset is removed 
by taking the difference between the two time traces acquired on each sequence. 

The spin energy relaxation time Ti being ~ 0.4 s (see Fig.2e of the main text), we choose a 
repetition rate 'yrep sufficiently slow to allow full relaxation of the spins in-between successive 
sequences. For example, the spin-echo spectroscopy shown on Figure 2 is acquired with 
Irep = 0.04 Hz, the data of Figure 3b&c with = 0.1 Hz and the absorption data of 
Figure 3d with q^ep = 0.3 Hz. 

II. BISMUTH DONOR SPIN AND COUPLING TO THE RESONATOR 

Neutral bismuth donors in silicon have a S' = 1/2 electron spin and a nuclear spin 1 = 9/2 
that are strongly coupled by an isotropic hyperhne interaction term A/2t[ = 1.45 GHz. The 
system is described with the following Hamiltonian, [28], where qe/27r = 28 GHz/T and 
qn/27r = 7 MHz/T ; 

^/h = B-(qeS(g)l-q„l(g)I) + AS-I (SI) 

In the limit of low static magnetic held Hq, the 20 electro-nuclear energy states are well 
approximated by eigenstates of the total angular momentum F = S -|- I, which can be 
grouped in an F = 4 ground and an F = 5 excited multiplet separated by a frequency of 
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bA/2T[ = 7.35 GHz in zero-field, shown on Figure 1. For a given low static field Bq oriented 
along z, only transitions verifying \Amp\ = 1 have a sizeable matrix element (equals to 
the Sy matrix element) and so may be probed with an excitation held orientated along x 
(or equivalently along y). We give in the table below details on the two transitions that are 
accessible to our resonator. 


Transition 

Expected crossing field 

dfIdB 

1 (mp 1 1 m'p^ 1 = 1 (mp \ Sy \ m'p^ \ 

mp = —4 —)■ —5 

5.16 mT 

-25.1GHz/T 

0.47 

mp = —3 —)• —4 

6.68 mT 

-19.2GHz/T 

0.42 


TABLE II. 


Single-spin conpling to the resonator 

In the experiment, the static magnetic held Bq is applied parallel to the surface (see 
Fig. S2) along an axis Z that can be decomposed along the orientations dehned in Fig. S3 
as : 


Bo = HqZ = Ho cos(6*)z-|-Ho sin(6*)x. (S2) 

A full orthonormal basis is provided by the combination of X, Y, Z, with X = cos^x — 
sin 9 z, and Y = y. The total magnetic held B is the sum of the static bias magnetic held Bq 
and of the microwave held generated by the resonator Bi = (5B(ac + aJ), where we introduce 
the magnetic held rms huctuations at the spin location (5B and the resonator annihilation 
(resp. creation) operator he (resp. al). The held Bi is located in the plane perpendicular to 
the resonator wire; as a result one can write hB = hH^jX -|- SByj = SBx^ + SByY + SBz'Zi, 
with 6Bx = cos 96Bx, SBy = SBy and 6Bz = sinOSB^;. 

Projecting the total Hamiltonian of a single Bismuth donor in the held B on the Hilbert 
space spanned by the two levels Imp), \m'p) and introducing the usual Pauli operators yields 


H/h = ujQa\ac 



+ 7e [{mp\Sx5Bx + Sy5By\mp)a+ + {mp\Sx5Bx + Sy5By\mp)a-] (he -1- a\). 


(S3) 
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(In this equation we left aside the Sz terms as they contribute only as negligible 
fast-rotating terms). As noted in the previous paragraph, {mF\SxSBx + SySBYlrri'p) = 
{mF\Sx\m'p){SBx + iSBy), where {mF\Sx\m'p) values are shown in Table II (~ 0.5). Af¬ 
ter performing the rotating-wave approximation, the full system Hamiltonian takes the 
Jaynes-Cummings form 


H/h = uoaldc —-|- 5 f(e*'^°(j+ac -|- e *'^°d_a]!), (S4) 

where 00 is an irrelevant phase that can be absorbed by a re-definition of the energy 
levels, and g is the spin-resonator coupling constant given by 

g = {mF\Sx\m'p)'ye\j5Bl + {cosdySBl. (S5) 

The vacuum field fluctuations JB have a spatial dependence, fixed by the shape of the LC 
resonator mode, which implies that the coupling constant to the resonator will also follow 
the same spatial dependence. We determine this spatial dependence numerically in the 
following way. First, the spatial distribution of the current fluctuations in the resonator wire 
is computed, knowing that the integrated current over the wire cross-sectional area is given 
by 5i = ujq^/KJ2^, Zq = \JL/C being the resonator impedance estimated to be AAVt using 
the electromagnetic simulator CST Microwave Studio. For our 50nm-thick aluminum films, 
the current density is assumed constant in the y direction with an x-dependent integrated 
value SJ{x) given by [32]: 


6J{x) = { 


(5J(0)[1 — (2x/wY] for |a;| < — A^/(26)| 

(5J(^tc) exp — [(Itc — |a;|) 6/A^] ioi \\w — \^/{2h)\ < \x\ < \w (S6) 
(1.165/A)(m;&)1/2^J(0) 


for X = \w. 


In these expressions, w = 5 pm is the width of the wire, 6 = 50 nm is its thickness and 
A = 90 nm is the penetration depth for our A1 film. The normalization constant (5J(0) 
is determined by the condition that J^^i^5J{x)dx = 5i. From the current distribution, 
the spatial dependence of JB is readily obtained using Comsol, and is shown in Fig.S3. 
Importantly, we note that the Bi field is essentially along x in the region just below the 
wire, and essentially along y in the region immediately outside of the wire. According to 
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FIG. S3, (a) Spatial distribution of the current rms vacuum fluctuations flowing through the 

resonator inductance, corresponding to an impedance of 4411. (b) Scheme of the resonator, with 
corresponding directions used in the text, (c)rms vacuum fluctuations of the magnetic field at the 
given red cross-section on b (d) x (red) and y (blue) components for the vacuum fluctuations of the 
magnetic field at y = —lOOnm (left axis) and for the microwave field Bi corresponding to an input 
power Pin = —88.5dBm = 1.4pW (right axis), (e) Spin-echo spectroscopy realized for Bq = Bq-z 
(blue circles) and Bq = Bq • x (green circles) allowing to make the distinction between spins lying 
next (strong SBy) and under (strong 5Bx) the aluminium wire. 


Eq.(5), one thus expects the coupling to the microwave held to be strongly ^-dependent for 
donors located immediately below the wire, and to depend negligibly on 9 for those outside 
of the wire. 
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Sub-structures of the lines profiles 


The spectroscopy of the transition mp = —5 —)■ mp = —4 was realized with the magnetic 
held Bq applied snccessively along two directions : along z (6^ = 0°, parallel to the wire), 
and along x (6* = 90°, perpendicnlar to the wire). As shown in Snpplementary Fignre S3, 
the low-held snb-peak appears only for 0 = 0° whereas the npper-held strnctnre remains 
nnchanged for both orientations. 

The prohle of the microwave excitation held components and Biy at a depth of 
100 nm (implantation prohle peak) is shown in Fig. S3d. As mentioned earlier the held 
below the wire is essentially along x, whereas it is essentially along y ontside the wire. 
When Bo = Bqz, Bi is transverse to Bq for all spins. When Bq = 5ox, Bi is transverse 
to Bq only for spins ontside the wire. This strongly snggests that when Bq = Bqz both 
spin families can contribnte to the signal and when Bq = Sqx only spins ontside the wire 
contribnte. 

Al wire e X 10-4 



FIG. S4. 

As a resnlt of the larger thermal expansion coefficient for alnmininm than for silicon, 
once cooled the region nnderneath the wire experiences a diherent strain to the region im¬ 
mediately ontside of the wire see Snpplementary Fignre S4. Silicon is an indirect bandgap 
semiconductor with a six-fold degenerate condnction band minimnm. The sharp conhn- 
ing potential of a donor canses the six-fold degenerate gronnd state to split into a sin¬ 
glet gronnd state: Ai = ;^{1,1,1,1,1,1}, a triply-degenerate excited state: T 2 (x,y,z) = 
;^{1, —1, 0, 0, 0,0}, ;^{0, 0, 1 , —1, 0, 0}, ;^{0, 0, 0, 0, 1 , —1} and a donbly-degenerate excited 
state: F^(xy,xyz) = ^{1,1, —1, —1, 0, 0}, “1; 2, 2}. Strain has the effect of low- 










ering the energy of conduction band minima (or “valleys”) in the direction of compressive 
strain and raising the energy of valleys in the direction of tensile strain. Each state will 
therefore have a shift in energy that depends on the applied strain and its valley composi¬ 
tion [33]. In addition, strain mixes the ground state Ai with the doublet excited states E^y 
and Exyz. 

The I = 9/2 nuclear spin of bismuth means that it possesses an electric quadrupole 
moment Q. The quadrupole moment can interact with an electric held gradient (EFG), 
which, for example, can be generated by the electron wavefunction through the operator [34]: 


(S7) 

rrEFG _ 6 3 q;/3 — 

«/3 - ^5 

where a oi 13 = x,y,z are the crystal or principal coordinate system. In the absence of 
strain, the ground state electron wavefunction is the Ai state, which is a symmetric com¬ 
bination of the six valleys. This symmetric state produces no EFG and thus a vanishing 
quadrupole interaction (QI). On the other hand, the excited states are an asymmetric com¬ 
bination of valleys and result in an asymmetric charge distribution and a non-zero EFG. For 
the case of strain applied along the z principle axis, an EFG is produced through mixing 
with the doublet excited state E^y^. The quadrupole coupling for a held Bq applied in the 
direction of this EFG is given by the interaction Hamiltonian: 

^Qi = eQKz ( 3 ? -?(? + !))/ (4/(2/ - 1)) (S9) 

From Equation (S9), it is evident that the QI produces an energy shift only on transitions 
whose states have differing mi. This is applicable to our low magnetic held transitions, which 
are highly mixed in the electron-nuclear basis. The sign of the EFG - and consequently the 
quadrupole shift of the spin transitions - depends on the sign of the induced strain, which 
as shown in Figure S4 is opposite for donors underneath the wire and to its side. 

The measured spin resonance lines in our device are split above and below the theoretical 
held values (see Figure 2 of main manuscript), suggesting an underlying interaction with 
both positive and negative frequency components. As we have shown, such a frequency 
distribution could be explained by a strain-induced QI. This is completely consistent with 
the observation above that the two sub-peaks are subject to Bi helds of diherent orientation. 
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FIG. S5. (a) Amplitude JPA signal-to-noise ratio SNR = ^Pout/Pnoise as a function of y/G for an 
input microwave tone, measured with a Spectrum Analyzer, showing a six-fold improvement with 
the JPA on. The JPA is operated in the non-degenerate mode. Open blue circles are experimental 
data, red curve is a fit as explained in the text, (b) Echo signal-to-noise ratio (in amplitude) as a 
function of y/G (blue open circles), with error bars estimated as explained below. Red dashed line 
shows the amplifier SNR{y/G) curve, rescaled to these data, showing that the spin-echo signal-to- 
noise is increased as expected from the JPA, until saturation. 

Spins constructing the high field peak experience a Bi that is always perpendicular to Bo, 
which is true of spins to the side of the wire. Due to strain these spins see a different QI, 
which would explain such a shift in the resonance held. A quantitative comparison between 
theory and experiment is the subject of ongoing work. 


III. EXPERIMENTAL DETERMINATION OF THE SIGNAL TO NOISE RATIO 

We now discuss the experimental determination of the SNR. We hrst characterize the 
signal-to-noise ratio improvement brought by the JPA itself over the following HEMT am- 
pliher cooled at 4 K. A continuous microwave signal at frequency u is sent directly on the 
JPA and its output spectrum Pout is measured with a spectrum analyzer for various JPA 
pump power settings in the non-degenerate mode. The JPA gain in the different setting may 
then be computed by G = Pout/Pout off). The same experiment is then repeated with¬ 
out any input signal so as to obtain Pnoise^ the noise power in the measurement bandwidth 
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of 100 kHz. The amplitude signal-to-noise ratio is then evaluated by SNR = \/Pout/Pnoise 
and is shown in Fig. S5a. 

It follows the expected dependence SNR oc a/G /{{G — l)n + Usyst)-, n being the total 
noise photon number with the JPA ON as dehned in the main text, and Usyst being the 
number of noise photons added when the JPA is off. This yields a ratio Usyst/n = 36 as 
observed in similar setups [20], indicating Usyst ~ 36 and Uamp+neq ~ 1 therefore approaching 
the quantum limit, with equal contributions from Ueq and Uamp (while Ug — O.Oln has a 
negligible contribution as verihed experimentally). 

We then study the SNR of the spin-echo at Bq = 51.8 Gs for various JPA gains, using 
homodyne demodulation. For that we first choose the local oscillator phase such that the 
echo is entirely on one quadrature (J). We then average 10 spin-echo signals yielding time- 
traces I{t) from which we compute the integrated echo amplitude S = I{t)dt-, 

^ P{t)dt when the microwave pulses are off. 


the noise is then obtained as Af = m j, ^ jp 
The noise being determined with an averaging of 500 traces, the statistical uncertainty on 
the SNR comes from the signal, so that the absolute uncertainty esNR = ^s/Af = 1/vTo 
since signal traces were averaged 10 times. As shown in Fig. S5, the resulting SNR = S/Af 
follows the same dependence as already obtained for the JPA itself, which shows that with 
the JPA on, the detected spin-echo reaches the quantum limit of sensitivity, with a gain in 
sensitivity by a factor x7 compared to JPA OFF. 

It is possible to further improve the SNR by using the JPA in its degenerate mode, by 
setting Up = 2uo and setting the relative phase of the pump tone such that the amplihed 
quadrature is I. This increases the gain by 6 dB while only increasing the noise power by 
3dB. The expected increase of SNR by \/2 is indeed approximately observed, yielding a 
absolute signal-to-noise ratio of 7, a factor xll larger than with the JPA OFF as shown in 
Fig.3 of the main text. 


IV. NUMERICAL SIMULATIONS 

The goal of this section is to detail the numerical hts of the absorption and the spin-echo 
sequences shown in hgure 3 of the main text. From this simulation, we extract the number 
of spins probed in a single spin-echo sequence as well as the measured spin concentration. 
In order to reproduce the spin dynamics, the inhomogeneity in both spin frequency and 
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coupling strength is taken into account by dividing the ensemble into a sufficiently large set 
of homogeneous sub-ensembles and integrating the equations of motion for the resonator 
held and the spin components of all of the sub-ensembles using the model already described 
in [35-37], 

Model 

Consider an ensemble of N spin-1/2 particles of frequency Uj. Each spin couples to 
the resonator held (described by creation and annihilation operators aj and he) with a 
coupling constant Qj and a Jaynes-Cummings interaction (see Eq.4 above). The total system 
Hamiltonian is then 


^ N N 

H/h = uoalcLc + - + i\/ Ki/2(/3aj] - fi*a^ -F (SIO) 

i=i i=i 

with the Pauli operators of spin j for k = , 2 :}, and [3 the amplitude of the 

microwave held driving the cavity input in the laboratory frame. The equations of motion 
are then integrated under the Markov approximation to incorporate the ehect of resonator 
leakage and spin decoherence [37]. This numerical simulation yields the dynamical evolution 
of the mean values of the resonator held quadratures as well as of the spin operators. 

The inhomogeneity in both spin resonance frequencies and coupling strengths is taken 
into account by dividing the entire inhomogeneous ensemble into M homogeneous sub¬ 
ensembles, Ml, M 2 , ■ ■ ■ ,Mm, each of them describing spins having an identical frequency 
Um and coupling to the cavity held gm- For a sub-ensemble m, we dehne the total number of 
spins as Nm and the three spin collective operators as with i G {x,y, z}. 

Spin decoherence is treated by including a spin dephasing rate 7 _l = I/T 2 and a spin 
energy decay rate yy = 1/Ti. We use the experimentally measured coherence time T 2 = 9 ms 
(see Fig.2 of the main text). The energy relaxation rate on the other hand is dominated 
by Purcell relaxation through the cavity, meaning that Ti is longer for spins detuned from 
the cavity than for spins perfectly at resonance as will be discussed in later work. This is 
captured by dehning for each ensemble as k , with = oom — ojq. Note that the 

relaxation time shown in Figure 2 Ti = 0.4 s was taken with a very narrow bandwidth pulse 
so as to obtain only the contribution from the spins that are on resonance with the cavity. 
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This leads us to introduce for each sub-ensemble an effective initial polarisation {t = 
0). Indeed, every experimental sequence is repeated several times at rate {irep ~ 3 to 
10s) and the results are then averaged. This waiting time is long enough compared to 
Ti to be neglected for spins at resonance, however detuned spins have a longer Ti and thus 
do not fully relax between two consecutive sequences, contributing less to the signal than 
spins at resonance. To take into account this effect, we define an effective initial polarisation 
= 0) for a sub-ensemble depending on its relaxation time: 

= 0) = -N^ X (1 - (Sll) 

The first step before performing the simulations is to determine in the context of the 
experiment the size Nm of each sub-ensemble, which requires knowledge of the distribution 
of coupling constants and resonance frequency within the spin-ensemble. 

Determining the coupling constant distribution 

From the simulation of the vacuum fluctuations of the magnetic field shown in Figure 
S3, one can compute the coupling constant distribution using Eq. S5. The measurements 
which we want to simulate were performed with the magnetic field Bq aligned along the wire 
in the z direction, on the low-field peak of the structure shown on Figure 2, at 5.18 mT, 
attributed to spins lying under the wire. As a consequence, we compute the distribution only 
for this subset of spins, imposing |a;| < 2.5/im and 0 = 0 in the formula. The inhomogeneous 
implantation profile of the spins is taken into account by appropriately weighing the coupling 
constant distribution. In order to normalize the distribution we take ^ p(fi'm) = 1- As shown 
in Fig. S6a, this yields a very asymmetric distribution sharply peaked around g/27r = 56 Hz. 

Determining the spin frequency distribution 

The cavity linewidth (20 kHz) is two orders of magnitude smaller than the spin linewidth 
3.25 MHz. In order to avoid numerical errors, we use a sampling of 1 bin per IkHz, over a 
range of 450 kHz. The zero order approximation would be to assume a square spin frequency 
distribution, nevertheless we introduce a tilted square distribution to take into account more 
precisely the shape of the line. Figure S6. The relative slope is derived from the observed 
spin-echo signal. At 51.8 mT, this yields a tilt of 10% on the chosen range. 
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FIG. S6. (a) Coupling strength distribution extracted from magnetic field simulation for |x| < 
2.5 fim, with Zq = 440, weighted by spin concentration and normalized to unity (green line). The 
black circles show the discrete distribution used in the simulation, with Mg = 50. (b) Tilted square 
distribution used in the simulation, with Miyi^s = 450 (black line). The resonator transmission is 
plotted for comparison (blue line) 

Evaluating the number of spins 

In order to determine the absolute scaling of the spin distributions, we measure the time- 
dependent absorption of the spins, as shown in Fig. S7. A hrst 500 ps-long pulse P of 
power Pin is sent to the unsaturated spins, leading to the absorption of the signal (Fig. S7c) 
whereas a second pulse of same power Pin sent immediately after a strong microwave pulse 
whose role is to saturate the spins shows only the cavity dynamics (subset b). The sequence 
is repeated 1000 times with a repetition time = 3 s= lOTi to let the spins relax back to 
their ground state in-between the experimental sequences. 

The transmitted pulse P shows two prominent features : the Rabi oscillation transients 
at the beginning, which are characterized by an oscillation frequency Qr, a decay time and 
an initial amplitude, and the free-induction decay (FID) of the spins which gives rise to the 
emission of a microwave signal even after the cavity held has decayed. 

To simulate such a sequence, the coupling distribution and the spin frequency distribution 
detailed above have been scaled with a total number of spins Ntot = The cavity 

m 

parameters ki, K 2 , kr, and Uc as well as spin relaxation rates 7 _l and yy are experimentally 
determined as explained earlier. The simulation of the saturated pulses transmission is in 
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FIG. S7. (a) Absorption sequence consisting of a first 500-us-long pulse at power P(nphotons), 
followed by a strong saturating microwave pulse immediately followed by a second 500-us-long 


pulse at same power P(nphotons). (b) Saturated pulses taken with average number of photons Uph 
for the intra-cavity field, rescaled to same amplitude with an additional offset and averaged 1000 
times. Open circles : data, solid lines : ht. (c) Absorbed pulses, for same Uph, rescaled thanks to 


saturated curves and averaged 1000 times. Open circles : experimental data, solid lines: fit 

quantitative agreement with the data, without adjustable parameter (Figure S7b). 

The simulated absorbed pulses are calculated with the same cavity parameters and an 
overall distribution scaling factor Ntot = 2 x 10^. This parameter gives a good agreement 
with the FID part of the signal as well as with the Rabi oscillations amplitude. To obtain the 
right Rabi oscillation frequency, we simply need to scale by a factor 77 = 1.1 the power at the 
cavity input, compared to the input power estimated from the cables and filters attenuation. 
With only those two adjustable parameters, we are able to reproduce quantitatively the spin 
absorption (Figure S7c), which is a hrst validation for using this model as an evaluation of 
the number of spins contributing to the signal. 

This approach neglects one aspect of the experiment. In our sample, the inhomogeneous 
broadening of the line is caused by strain due to the presence of the aluminium wire. One 
can thus expect a correlation between the frequency of a given spin and its coupling to 
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the resonator. Taking this correlation between coupling constant and spin frequency into 
account quantitatively would however require a microscopic modelling for the strain which 
is not available so fair. The excellent overall agreement between theory and measurements 
indicates that the error is minor, although this approximation may account for the small 
discrepancy between simulated and experimental Rabi oscillations that can be noted for low 
powers in Figure 4c. 


Evaluating the number of spins contributing to the spin-echo signal 

Having as explained above calibrated the absolute scale of the spin distribution, we 
now evaluate the number of spins involved in a spin-echo by simulating a full Hahn echo 
sequence, keeping exactly the same parameters for the spin distributions. The input power 
of the simulated 7r/2 and tt pulses are calibrated by simulating Rabi oscillations. We hnd 
that in the simulation the tt pulse power is only 1 dB away from the experimental one, which 
further confirms the validity of our model. 

The spin echo sequence was acquired with the JPA off, in order to avoid its saturation by 
the drive pulses which would distort them. The output amplitude is scaled by comparing the 
theoretical and experimental decay of the two excitation pulses; with only this adjustment 
factor the simulated echo is found to be in quantitative agreement with the experimental 
data as shown in Fig.3 of the main text. 

To evaluate the number of spins excited during the spin-echo sequence, we extract from 
the simulation the time-dependent mean spin polarization {Sz), as shown in Fig. S8a. We 
consider more particularly that the quantity {Sz(t > n/2)) — {Sz(t = 0)) is a direct estimate 
of the number of spins excited by an Hahn-echo sequence. After the exponential decay of 
the 7r/2 excitation pulse, this value increases by 1.2 x 10^. We thus come to the conclusion 
that 1.2 X 10^ spins participate to the echo shown in Fig.3 of the main text, which is detected 
with a SNR = 7, which yields the sensitivity reported in the main text. 


Evaluating the concentration of bismnth donor spins 

Thanks to the determination of the spin density in absolute scale, we can compare the 
estimated number of spins to the known number of implanted atoms, to check that the two 
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FIG. S8. (a) Red: Simulated 7r/2 excitation pulse. Blue: {Sz{t > vr/2)) — {Sz{t = 0)), allowing to 
extract 1.2 x 10'^ excited spins during this excitation pulse, (b) Spin density profile (red circles) 
oc X ^TT (defined in the main text) rescaled to f p{ijjspins)dojspins = 1 (grey filled area). 

are consistent. We can scale the spin density p{(jJs) oc Ag x thanks to the absorption 
simnlation. By integrating the lower field peak, Fignre S8b, we find a total nnmber of spins 
of 1.07 X 10® contribnting to the absorption. As Bismnth donors have a nnclear spin of 
9/2, this represents only one-tenth of the total amonnt of spins. Also as this peak has been 
identified as signal emitted only by spins nnder the wire, the amonnt of spins is dilnted in a 
snrface S'^jre = L-wire X W^ire = 720pm X 5pm = 3.6 x 10“® cm^, leading to an experimental 
snrface concentration [Bi]Exp = 10 x 1.07 x 10®/(3.6 x 10“5) = 2.97 x 10^^ cm“^, a nnmber 
which is only a factor 3 lower than the snrface concentration extracted by SIMS measnrement 
[Bi]siMS = 9-45 X 10^^ . 

This ratio of 30% can be explained by two factors: hrst that the activation of bismnth 
atoms i.e. the migration of bismnth atoms from interstitial implantation site to snbstitn- 
tional site by rapid annealing is not total. This factor has been evalnated in [29] to be 60% 
by electrical measnrement. The additional factor 2 with onr experiment may be dne to a 
fraction of bismnth atoms being in an ionized state and thns not contribnting to the ESR 
signal. 


Approximate analytical formula for the expected spin-echo signal-to-noise ratio 

We now derive an approximate analytical formnla for the pnlsed ESR spectrometer sen¬ 
sitivity, which is given in the main text. The goal is to provide an analytical estimate for 
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the experimentally obtained sensitivity and thus identify how the sensitivity scales with the 
physical parameters. 

Compared to the model used in section IV, we consider here for simplicity that all spins 
have equal coupling constant g to the resonator. As a result in the frame rotating at the 
resonator frequency uo, the Hamiltonian is: 






(S12) 


3 " 

To obtain analytical results, it is convenient to assume that the spin detunings Aj are 
distributed according to a Lorentzian function /(A) = with a FWHM of w. More 

realistic distributions would only change the hnal results by a factor of order unity. For 
simplicity again, we neglect decoherence of the spins. The following equations of motion 
then describe the evolution of mean values: 


= (S13) 

3 

(S14) 

We dehne S- = \ and formally integrate (S13), 

{dc{t)) = —ig f l^S-{t'))dt'. (S15) 

J —OO 

From this equation we see that (he) attains at most a value of ~ gN/n, and if g times 
this value is much smaller than typical detunings, A^, he., the cooperativity parameter 
C = 1, the second term of Eq. (S14) can be neglected. We assume this is the case 

and we thus treat the spins as evolving freely between the pulses. 

We will focus on the case of a standard Hahn-echo sequence, n/2-T-n-T, and we assume 
perfect initial 7r/2 and refocusing n pulses at times —2r and —r, respectively. Prior to 
t = —2r, the spins are polarized with polarization p < 1 along the —; 2 -direction such that 
= —p and {dx'^) = {dy ^) = 0. The case p = 1 corresponds to a perfectly polarized 
sample with all spins in the ground state. When the spin ensemble is subjected to the 
perfect 7r/2 pulse around the x-axis at time t = —2r, a state with (dy'^^(O)) = p and hence 
{d^l\o)} = — y is prepared. The refocusing tt pulse at t = —r implies that a spin echo will 
occur at time t = 0, where the free evolution, {d_ (t)) = — —, cause all spins to be in 
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phase. The resulting total spin (S'_(t)) can be calculated in the continuum limit, where the 
sum over spins is replaced by an integral over detunings, 

{S-(()) = J f(A)=^e-*^‘dA = (S16) 

Inserting this into Eq. (S15) leads to the following time dependence of the resonator held: 


(flcW) 
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t < 0, 
t > 0. 


(S17) 


Now, as indicated in Fig. (1) of the main text, this resonator held is coupled with rate K 2 
into a Josephson parametric ampliher (JPA). The input held, Sin, to this ampliher is related 
to the resonator held. Sc, by (Sin(t)) = yjl^{ac{t)) and is normalized such that (S|j^(t)Sin(t)) 
represents the number of microwave photons per time incident on the ampliher. We shall also 
dehne the input held quadrature variables by Xin(t) = and yin(t) = ^ 

Since (Sc(t)) is real-valued in our calculations, see Eq. (S17), the mean signal is carried by 
the X-quadrature. 

So far we have only considered mean values of the spins and the cavity held, and 
we shall apply an operator description of the amplihcation stage to assess the noise on 
the measurement signal. To this end we dehne single modes of the propagating held by 
din = / ain(t)u(t)dt, where u(t) is a mode function chosen to be real valued and fulhlling 
J[u(t)]‘^dt = 1. In this case we have [ain,h|^] = // u(t)u(t')[din(t), dl^(t')]dtdt' = 1 due to the 
free held commutator relations [ain(t),ain(t')] = — t'). The corresponding single mode 

quadrature variables, Wn = / Xi^{t)u{t)dt and bin = ^ %nif)u{t)dt, fulhll [Xin,I1n] = 
and the minimum uncertainty state, obtained for the vacuum state or coherent states of the 
held, must obey (AXj^) = = f. 

Noise in the propagating held is conveniently characterized by its dimensionless power 
spectrum S{u)/hu = (AAj^ -|- AY^^^. At the cavity output, noise arises from both electro¬ 
magnetic equilibrium huctuations, characterized by Ueq, and from possible extra noise due 
to spontaneous emission of the spins, as will be discussed further, with a contribution Usp. 
Our experiments being performed at temperatures such that kT huo, the electromagnetic 
held at equilibrium is indeed very close to its ground state so that we can safely assume that 
Ueq = 1/2. In total we get (AA/n) = l(nsp -1- 1). 

When the signal pulse, emitted from the resonator, is transmitted through the ampliher, 
its power is increased by the gain G. However the amplihcation process itself can add extra 
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noise to the output field, characterized by a dimensionless power density namp, and further 
degrade the signal-to-noise ratio. Following Caves [19], two cases should then be envisioned 
to describe the statistics of the field at the amplifier output. If the amplifier is in the so- 
called non-degenerate mode, its single mode output is described by the field annihilation 
operator: 

hout = din -|- VG — 1 b\^. (S18) 

To ensure the bosonic commutator relation of the amplified signal operators, an idler mode 
operator, obeying [bid, ~ must be included in the amplifier relation [19]. It should also 
be noted that the gain G is generally a function of frequency, which may distort the temporal 
shape of an amplified pulse. However, the measured 1.7 MHz wide gain profile, see Fig. 1(f) 
of the main text, supports that we can assume a constant gain across the bandwidth of the 
signal and hence for the single mode defined by u{t). 

We introduce quadrature operators for the output and idler modes in a similar manner 
as for the input field, and the mean value of the input field is simply amplified as (Xout) = 
'/G{Xia)- However, the noise in the output has contributions from both the signal and the 
idler mode, 

(AXL) = + (G - l)(AXf,). (S19) 

Assuming the idler mode thermalized at a temperature T yields (AX;^) = |(1 -|- 2h), n = 
1/(e"^ — 1) being the mean thermal photon number. At high temperatures 3> 1 so that 
n ~ the thermal state of the idler yields an overwhelming contribution (G — 1)^?:^ to 
the readout noise on the X quadrature, with rZamp = h- This is in particular the case in our 
experiments when the JPA is turned off and the signal is exclusively amplified by the HEMT 
amplifier at the 4K stage, for which namp ~ 50. When the JPA is on, this contribution is 
minimized since amplification is now carried on at a temperature T such that k^T hujQ, 
with an amplifier that reaches the quantum limit. In these conditions, namp = 1/2. 

In total, we find that the amplification obeys the following relation for the signal-to-noise 

(Wut)7AW ^ G(n,, + ,v) 

(Xm)^/AXj^ G(neq J- ngp) -|- (G — l)namp 

In the special case that namp = ''T-eq = \ and assuming ngp = 0, this equals which 

yields, in the limit of large G, the well known factor of two (i.e. 3 dB) reduction in squared 
signal-to-noise by phase insensitive amplification with no excess noise. After amplification 
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by the JPA, the signal level is sufficient that amplihcation and homodyne demodulation do 
not further degrade the signal-to-noise. In our analytical estimate of the ESR sensitivity, 
we shall proceed with the assumption that the temperature is sufficiently low, the gain 
is sufficiently high, and the excess spin noise is negligible (we shall return briefly to an 
assessment of this assumption). 

The other case to consider is the one where the ampliher is phase-sensitive, with 
quadrature-dependent gains Gx,y- This ensures the correct output commutator relations 
without requiring the addition of idler mode noise as explained in [19]. An ideal ampliher 
at the quantum limit may then verify Gx = Gy^ 3> 1 while implementing noiseless am- 
plihcation of one of the two quadratures, implying that in this case = n/2, with 

^ '^eq,X T ^amp,X &Ild UampjX 0- 

Summing up the discussion, one sees that the noise on quadrature X referred to the JPA 
input can be written as = n/2, with n = neq -|- Uamp + ^^sp- In our experiment, 

Ueq ~ 1/2, whereas Uamp ~ 50 if the JPA is off, Uamp ~ 1/2 if it is operated in the non¬ 
degenerate mode, and Uamp ~ 0 if it is operated in degenerate mode. 

Eqs. (S19, S20) account for the noise and signal-to-noise for measurements of the con¬ 
tinuous output amplitude signal weighted with u(t). The optimal choice for u(t) is the 
one that maximizes the weighted signal (Xout) without altering the noise (assuming uni¬ 
form broad band noise of the idler mode bid)- In an experiment one may choose u{t) 
as the measured shape of the emitted pulse averaged over many experimental runs. As 
shown by the red curve in the inset of Fig. 3(c) in the main text, we can obtain the 
same shape by a numerical calculation of the spin dynamics. Since we are here inter¬ 
ested in analytical estimates, we choose the result given by u{t) = (Xin(t))/(Xin), where 
{Xin(t)) = |(ain(t) + a[jj(t)) = ^^{cLcit) + alif)) = ^/K^{ac{t)), with an explicit expression 
given in Eq.(S17). To obtain the correct normalization J[u{t)]‘^dt = 1, we calculate the 
squared signal as: 


(Ain)' 


{xut))^dt 


2 g‘^p^N‘^K 2 {K -I- 2w) 
{k, -|- w)‘^wn 


(S21) 


The output signal-to-noise ratio then reads: 


|(Ao„t)| _ 2gpN lY I K2 {h + 2iu) 
AXout K, + w\n\ WK 


(S22) 


The minimum number of detectable spins, which dehnes the ESR spectrometer sensitivity. 
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is thus 


N ■ = 

-i ’ min 


K + W 


nwK 


‘^gp V 1^2(1^+ ‘^w) ‘^gp^ 1^2 ' (S23) 

The arrow indicates the limit of conventional ESR operation, w. Introducing the echo 
duration Te = w~^ as seen from Eq. S17, and taking the case of a critically coupled resonator 
for which n = 2 k, 2 , one obtains Nmin = which is the formula found in the main text 

in the experimentally relevant case p = 1. For our parameters, this yields Nmin = 400. 

This estimate can be further rehned to take into acconnt the fact that the actual width 
r/27r ~ 1 MHz of the spin frequency distribution is much larger than k. In this case, the 
7r/2 and tt pulses excite a subset of spins, and according to the numerical simulations, this 
subset has a Lorentzian prohle with w/2'k ~ 25 kHz, which is very close to k/2t: ~ 22 kHz. 
Hence, in the above expression (to the left of the arrow) one can take k ^ w ^ 2 k2 for 
critical coupling, which yields Nrr,\„ ~ ~ 3-10^, using g/2Ti ^ 55 Hz and p ^ 1. 

This rehnement should however not be considered too seriously when compared to the 
experimental data, given the other approximations that were made, snch as perfect tt pulse, 
Lorentzian line prohle, and optimal weighing fnnction u{t). 

The signal is farther enhanced, with no increase in the noise, when accumnlated over the 
CPMG echoes, i.e., by choosing the corresponding multi-peaked u{t), cf. Fig. 4(a) of the 
main text. In the analysis of the experiments, we assumed a weighting of the signals by 
tophat pulses of equal weight, and we, indeed, observed an improved signal-to-noise over 
the single pulse analysis. Let us here estimate the theoretical limitations of the CPMG echo 
spectroscopy, assuming a gradual, exponential redaction of the echo amplitudes. With m 
pulses in total, we may dehne the corresponding normalized mode function implicitly 

as: 


nw 


u 




(S24) 


(W) 

where is the outpnt quadrature for a single tt pulse and we assume that the echoes 

are non-overlapping and simply repeated with a period of T but damped by the rate T^pmg- 
This Tcpmg Kiay include any experimentally determined damping effects, snch as non-ideal 
vr pnlses causing a degradation of the signal. As in the above analysis, the normalization, 
J[u^^\t)]'^dt = 1 is automatically accounted for when calculating the squared signal as: 

r 1 fi—2.raT/T cpmg 

XPf = / (S25) 


i=o 
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The use of multiple pulses in the CPMG protocol only makes sense if T Tcpmg-, and 
the denominator in the above expression can be approximated as 2T/Tcpmg- Hence, the 
signal-to-noise ratio, SNE^ for m pulses becomes: 


SNR„ 

SNRi 


T 


CPMG 

2T 


— (^-'2'mT/TcpMG'j^ 


(S26) 


which behaves as y/m for a small number of pulses with 2mT Tqpmg and saturates at 
\/Tcpmg/‘^T when 2mT 3> Tcpmg- 

Let us hnally return to our assumption that the spin noise is negligible. To assess this 
issue we have performed a numerical calculation similar to the one shown in Fig. 3(c) in 
the main text, but including also the quantum noise using the methods of [35]. For an 
effective number of spins N ^ 1.2 ■ 10^, this calculation, indeed, shows an excess noise of 
~ 30 %, and we hnd that this noise is proportional to the number of spins involved. Since 
the excess noise is calculated for the situation where the signal-to-noise ratio is 7± 1, it must 
be seven times smaller, i.e., only at a few percent level when the signal-to-noise is unity. For 
this reason, it does not affect the experimental assessment of iVmin ~ 1-7 ■ 10^ in the main 
text. We also note that with the effective number of spins N ^ 1.2 ■ 10^, the cooperativity 
parameter reaches the value of C* = ~ 0.26. This number is also proportional to N and 

thus approximately seven times smaller in the case of a signal-to-noise level of unity, thus 
validating our assumption, C* 1, in the analytical estimate. 
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